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1 Introduction 



The framework introduced by Randall and Sundrum in ref. [1] has been playing an impor- 
tant role in high-energy physics in the last years. The Randall-Sundrum setup involves 
5D spacetime with the line element 

ds^ — a^{x^)dXfj_dx'^ — dx\ , (1) 

and the fifth dimension truncated at X5 = by the ultraviolet (UV) brane and ai — L by 
the infrared (IR) brane. The scale (warp) factor multiplying the 4D Minkowski background 
metric varies along the fifth dimension, which generates a hierarchy of scales of order 
a(0)/a(L) between the UV and the IR branes. For example, for the AdSs spacetime the 
warp factor varies exponentially and a huge hierarchy can easily be generated. 

The original motivation was to explain in this way the hierarchy between the Planck 
scale and the electroweak scale. However it has become clear that the scope of application 
is much wider. In particular, it is interesting to consider situations where the gravitational 
degrees of freedom can be decoupled. In this case one deals with 5D theories of gauge 
and matter fields in a fixed background of the form of eq. (1). This approach has turned 
out to be fruitful for constructing realistic models of the Higgs sector [2] , Higgsless elec- 
troweak breaking [3] and supersymmetry breaking [4]. Furthermore, AdS/CFT [5] applied 
to the Randall-Sundrum background [6] suggests that these models are dual descriptions 
of purely four-dimensional strongly coupled physics. There is also a connection between 
field-theoretical 5D gauge theories in warped backgrounds and the low-energy physics of 
pions and vector resonances [7], known as AdS/QCD. 

Gauge theories in D > A dimensions possess an intrinsic cutoff scale. The gauge 
coupling gn has dimension [length] -^Z^"^; therefore at high energies scattering amplitudes 
grow as E"^/^"^, leading to strong coupling in the UV. Perturbative computations have to 
be cut off below the strong coupling scale. In some phenomenological applications strong 
coupling occurs, in fact, not far from the TeV scale. Thus it is often desirable to consider 
a UV completion of higher dimensional theories so as to understand possible cutoff effects. 
Typically, this UV completion is assumed to be some sort of string theory. 

Deconstruction [8] is another option to model cutoff effects in higher dimensional gauge 
theories. It is a four-dimensional framework that typically involves a product gauge group 
and a set of bifundamental Higgs fields (the hnks). With an appropriate choice of 
representations and vev's of the links, such a setup, at low energies, reproduces the spec- 
trum and interactions of a higher-dimensional gauge theory with the gauge group G. The 
matching holds up to a certain deviation scale A/j, related to the magnitude of the link 
vev's. This deviation scale is identified with the cutoff of the higher- dimensional theory. 
Above Ad, deconstruction provides a UV completion in terms of a purely four-dimensional, 
weakly coupled and, possibly, renormalizable gauge theory dynamics. 

Deconstruction of 5D gauge theories in the AdSs Randall-Sundrum background was 
considered in refs. [9-13]. The warped fifth dimension was represented by a chain of 
bifundamental links <I>j with the vev varying as ~ . This setup turned out to be 
useful for clarifying several issues concerning the evolution of gauge couplings in AdSs 
[10-12]. However, these studies suffer from certain limitations. While computations in 
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AdS5 can be performed (both at the tree level and at the loop level) using familiar Bessel 
functions, no analytical methods have been available so far to handle computations in 
the deconstructed model. In particular, the spectrum and interactions of the massive 
gauge bosons have been determined only numerically. Some analytical results have been 
obtained, but only for g <^ 1, in which case dcconstruction does not have an obvious 5D 
interpretation. For this reason the relation between dcconstruction and the AdSs theory 
was somewhat obscure. It was not even clear if dcconstruction could really reproduce the 
AdSs physics in its entire perturbativity range. Furthermore, because of these technical 
problems, the deconstructed AdS^ setup has not been really useful for phenomenological 
applications. 

In this paper wc fill this gap. We present analytical methods to handle deconstructed 
AdSs models. The tools are provided by the mathematical theory of q-difference equations 
and their solutions. One of its branches deals with the so-called q-Bessel functions, which 
generalize the ordinary Bessel functions. We show that the q-Bessel functions are appro- 
priate to describe the spectrum and interactions of the deconstructed AdSs models. With 
these methods at hand, the calculabilities of the AdS5 gauge theory and its deconstructed 
version stand on an equal footing. 

The technical results we obtain help to clarify the relation between the 5D warped 
gauge theories and deconstruction. It was argued in the previous works [10, 11] that de- 
construction realizes a position-dependent cutoff. With the new methods at our disposal 
we are able to make this notion more precise and specify the parameter range in which 
deconstruction adequately approximates the 5D theory. The position-dependent cutoff is 
reahzed in the following sense. The IR brane scattering amphtudes in the 5D theory are 
reproduced up to the scale Am = Aa{L). On the other hand, the UV brane t-channel 
amplitudes are reproduced up to Auy = Aa(0), larger by the factor a(0)/a(L). The non- 
trivial thing about the latter result is that the spectra of massive excitations in the two 
theories deviate at a much lower scale, of order Am. Matching of the UV brane t-channel 
amplitudes in the two theories holds in spite of the fact that the number and the couplings 
of the exchanged massive gauge bosons are different. Our results also imply that, up to 
the scale Auy, the holographic interpretation of the deconstructed AdS5 models is similar 
to that established in [6] for the 5D models. 

The paper is organized as follows. In Section 2 we review the necessary technical 
material concerning gauge theories in AdS5. In Section 3 we study the deconstructed model 
and derive analytical formulas for the spectrum and propagators of the gauge bosons. The 
detailed discussion of these results and their consequences is presented in Section 4. Section 
5 contains a summary and points at future applications of our results. Deconstruction of 
5D gauge theories with most general boundary conditions is discussed in Appendix A. 
Finally, Appendix B contains a detailed and self-contained review of the theory of q-Bessel 
functions. 
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2 Review of gauge theory in AdSs 



In this section we summarize basic properties of 5D gauge theories in the Randall-Sundrum- 
type background. The fifth dimension is an interval ( equivalent ly, the orbifold jZ-^) 
parametrized by G [0,L]. Two boundary branes are located at ^5 = (the UV brane) 
and at 2:5 = L (the IR brane). The gravitational background is that of the AdS5 geometry: 

ds^ — a^{x^T]fxvdx^j^dxi, — dxl 0(2:5) = e"*^^® , (2) 

where k is the curvature scale and i^, • • • = . . . 3 label the 4D coordinates. The presence 
of the warp factor generates a hierarchy of scales of order l/a^ between the UV and the IR 
brane, where = a{L) = e"'^^. In the original Randall-Sundrum model ~ 10~^^, so 
as to explain the hierarchy between the electroweak and the Planck scale. In our analysis 
we allow for arbitrary -C 1. 

In the following we will focus on the dynamics of gauge theories on such a background, 
while fluctuations of the 5D metric will be ignored. The practical reason is that decon- 
struction of gravity encounters certain technical problems [14,15], which we prefer to avoid 
here. Formally speaking, gravity can always be decoupled by taking the limit M5 00, 
where M5 is the 5D Planck scale. More precisely, gravity couples strongly at the IR brane 
at the scale Ag(L) ~ aiMpianck ~ o,LMf^ /k}/'^, while at the UV brane it becomes strongly 
coupled only at Ag(0) ~ Mpianck- We assume we always deal with energies far below the 
respective gravity strong-coupling scales, and concentrate exclusively on the dynamics of 
the gauge sector. 

We consider a 5D gauge fleld Am = {An, A^) propagating in the background of eq. (2). 
The quadratic part of the action takes the form 

Ss^ J d'x £ dx,^^ + ^^{d,An - d,A,y^ . (3) 

We choose the boundary conditions as 

95^^(0) = Oa^L) = A(0) = A,iL) = . (4) 

The Neumann boundary conditions for A^ allow the existence of a 4D zero mode, so 
that 4D gauge symmetry survives below the compactification scale. Alternatively, one can 
impose Dirichlet boundary conditions for Af^ (and Neumann for A5), in which case gauge 
symmetry is entirely broken at low energies. For completeness we review this case and its 
deconstructed version in Appendix A. 

2.1 Spectrum and propagators 

There are several ways of approaching higher dimensional theories. Most often, higher 
dimensional fields are traded for an infinite number of 4D Kaluza-Klein (KK) modes. 
Another approach relies on 5D propagators, which effectively sum up the contribution of 
the whole KK tower. Finally, it is sometimes useful to employ the so-called holographic 
approach, which consists in constructing an effective action for the boundary values of 
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the bulk fields. Below we review the necessary technical background for each approach. 
Of course, these are merely different methods of organizing computations, and physical 
conclusions must be the same, regardless of which approach is used. 



2.1.1 Kaluza— Klein approach 

The gauge field is decomposed into KK modes as follows 



n=0 



n=l 



rrir. 



(5) 



The modes are eaten by the massive gauge fields but it is not necessary to fix the gauge 
at this stage. The KK profiles /„ should be chosen such that the quadratic action (3) is 
diagonal in the KK basis: 



This is obtained when are solutions of the equation of motion 

dlfn - ^kd^fn + mla~^{x5)fn = , 
satisfy the Neumann boundary conditions 

and are normalized as g^^ — 1. The solution is given by [16] 
/o 



(6) 



(7) 



(8) 
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Ana'^{x5) 



k J \ka{x5) J 
with the normalization constant^ 

2 \ -V2 



\ k ) ^ \ ka(x^) J ^ \ k J ^ \ka{x5) J 



n > , (9) 



An — 



kaL ' 



-1/2 



(10) 



Thus, the 5D gauge theory is rewritten in terms of a tower of 4D vector fields 
with masses m„ (^Is's are eaten). The zero mode corresponds to mo = 0. Its profile 
/o is a constant, thus the massless mode couples to all matter with equal strength given 
by go = gb/y/L- For the massive modes, the mass spectrum is given by solutions of the 
equation 

Mxhfe) - ^» (x) <") 

^We thank K.y.Oda for discussions on this and subsequent formulas in this subsection. 
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which can be well approximated by 



m„ ~ 7rkaL{n — 1/4) . 



(12) 



The spectrum is approximately linearly spaced with the mass gap Mk^- ~ kai. The 
massive modes couple non-universally to matter, the coupling depending on both the KK 
number and the position in the extra dimension. In particular, the couplings to the UV 
and IR branes are determined by the boundary values of the KK profiles 



/n(0) 



2kg 5 



2kg5 



fn{L) 



(13) 



- 1 



- 1 



Approximating the Bcsscl functions we find |/„(L)p ^ 2kg'l, so that all the KK modes 
couple to the IR brane with approximately equal strength. On the other hand |/n(0)p ^ 
2kaLg1{T^k / 2mn) log~^(m„/A;) for m„ -C k and |/„(0)p ~ 2kaLg1 for ^ k. This shows 
that KK modes are localized toward the IR brane and couple very weakly to the UV brane. 
For this reason physics on the UV brane may remain perturbative at energies much larger 
than the KK mass gap. 



2.1.2 Position-space approach 

An alternative approach to quantum computations in higher dimensions relies on 5D prop- 
agators. Since 4D Poincarc invariance is preserved, it is convenient to work in a mixed 
representation of 4D momentum space and position space in the extra dimension which, 
for brevity, we call the position space representation: 

P^.Ap\ ^5, Z/s) = j d^xe'P-{TA^{x, x,)A,{0, y,)) . (14) 

At tree level the propagator is an inverse of the kinetic operator (Fourier transformed to 4D 
momentum space). We choose the gauge-fixing term Cgf — —^^(c?^^^ — ^95(0^(^5)^5))^. 
Then the propagator satisfies the equation 

A {i~P^ ~ 95{a'^i^5)d5)]Vfiu + (1 - ^/OPuP'^) PupiP^, ^5, Vb) = iVfipS{x5 - 1/5) (15) 

95 

and the Neumann boundary conditions dx^Pfj,„{p^, X5, 2/5)1x5=0 = d^^Puvip^, a^s, 2/5)U5=l = 0. 
The solution is given by [17] 

-PML/(p^ x^, y^) = {^^^ - P(p^ a;5, 2/5) + ^^P(pVC, 2:5, 2/5) (16) 

where 

mgl [^o(f)'^i(tofe)) ~ ■^o(fc)^i(fc^)][^o(fcfz)'^i(fc^) ~ '^o(fcfr)^i(Mk)^] 
2k a(^5)a(l/5)[ro(f)Jo(^) - ^o(|)Vo(^)] 

(17) 



P{P ,2:5,2/5) = 
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for < y5 (the propagator is of course symmetric under the exchange ^ y^). 
We will focus on brane-to-branc propagators defined as 

Pw(p') = P(P^ 0, 0) , P,r{p') ^ P{p\ L, L) , P^sip') ^ P{p\ 0, L) . (18) 

They read 



p yo{ 
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p \ ' 
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Pbb{p^) - v(P\ 1 ( p \ 7 (P\v ( P \ ■ (19) 



Consider the Euclidean momenta = —Pe- Pe ^ kai all propagators approximate 
those of a 4D massless gauge boson with the gauge coupling go — g^/y/Z: 

PM-pI) ~ Pir{-pI) ~ Pbb{-pI) ~ '4- ■ (20) 

Pe 

Above the KK mass gap the momentum dependence of the propagators crucially depends 
on the position in the extra dimension. For <^ pe <^ k one finds 



Pvv{-pI) 
Pmi-Pl) 



Pe M^) \og{2ke-ypE) p| 
aiPE ' 



,2 ' 



The propagator on the IR brane shows a fall-off. which is the same behaviour as for 5D 
flat spacetime. Propagation between the branes is exponentially suppressed, which is also 
a characteristic feature of higher dimensional theories above the compactification scale. 
But on the UV brane the fifth dimension is screened: the propagator exhibits the usual 
4D l/E"^ behaviour, up to a "classical running" encoded in the logarithmic form factor. 
Such behaviour persists all the way up to the curvature scale. At even higher energies, for 
Pe > k: 



Pe) ~ — /PK\ ~ — ) ^iR\ Pe) ~ „ „ ■, ^bb[-Pe) 



PeM^) pe' a,PE' ^'^^ ^V2^^^ ■ 

(22) 

At very high energies, larger than the curvature scale, the propagators have an energy 
dependence analogous to that of the propagators in the 5D flat spacetime. The UV and 
IR brane propagators differ only by the inverse warp factor multiplying the latter. 
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The connection between the KK and the 5D position space approach is given by the 
spectral formula 

m 2 \ ■ fny^bjfniyb) 

P{p ,x5,y5) = — o 5—- 23 

n "■ 

This shows that the position propagator describes collective propagation of all KK modes 
between and y^. The KK masses and wave functions arc given, respectively, by the 
poles and residua of the propagator. The position propagator is thus a very convenient 
object, which encodes information about the whole KK spectrum. Moreover, locality in 
the 5th dimension is explicit in this formalism. For this reason, it is particularly useful for 
the computation of scattering amphtudes of fields localized on 4D branes. 



2.1.3 Holographic approach 

In the holographic approach we single out the UV boundary value of the 5D field, treating 
it as a distinct variable from the bulk or the IR brane values. The latter degrees of 
freedom are integrated out, leaving a non-local effective action Scs for the UV boundary 
value. At the classical level, integrating out amounts to evaluating the 5D action on a field 
configuration that satisfies the bulk equations of motion and the IR boundary conditions, 
with the UV brane value left as the "low energy variable" . In the gauge ^ = 1, such 

a configuration can be written as 

^"('■-^'-W^-^-W. (24) 

where the propagator P(p^, 2:5, 7/5) was defined in eq. (17). Inserting this into the 5D action 
we obtain the effective action for An{p): 

1 f d^j) — — % 

We see that, at tree level, the self-energy of ^^^(p) is an inverse of the UV brane-to-brane 
propagator in (19). 

In principle, this procedure can be applied to arbitrary geometries, whenever the bound- 
ary value is for some reason a convenient low energy variable.^ However for the background 
of eq. (2) the AdS/CFT correspondence [5] applied to the Randall-Sundrum geometry [6] 
gives it a special meaning. The gauge theory we study here is dual to some 4D strongly 
coupled CFT, with a 4D gauge field weakly coupled to a conserved current of the CFT. 
Hence H^,^ represents the connected correlator of two conserved currents, and J^,, in 
the dual 4D theory. This allows us to understand the pecuhar features of the UV propa- 
gator described in the previous subsection. In particular, the form of P\jv in the regime 
ka^ <^ Pe k (see Eq. (21)) is fixed by conformal invariance of the 4D holographic theory. 

Cutting off AdSs with the UV brane is interpreted as an explicit breaking of conformal 
invariance in the 4D dual by a UV cutoff of order k. On the other hand, the presence of the 

^This is the case, for example, in models with the SM gauge fields living in the bulk and the SM matter 
localized on the UV brane [18]. 
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IR brane is interpreted as a spontaneous breaking of conformal invariance in the 4D theory. 
This results in a mass gap of order kaL and in a discrete spectrum of resonances that are 
identified (up to a mixing with ^4^) with the KK modes in the 5D picture. According to 
AdS/CFT, the tree-lcvcl approximation of the 5D theory corresponds to a large- A'^ limit 
in the dual theory. In fact, it is well known that, in the large- A^ limit, the exact two-point 
correlator function can be written as a sum of infinitely narrow resonances, in agreement 
with (23). 



2.2 Scales 

We close this section with a discussion of the energy range where perturbative gauge 
theories can be applied. As we are interested in the parametric dependence only, we do 
not display any numerical factors explicitly. 

So far we have encountered two scales: the curvature scale k and the KK scale Mkk ~ 
aik, that mark a qualitative change in the behaviour of the propagators. In 5D gauge 
theories the gauge coupling has dimension [length] therefore another scale l/gl appears. 
This quantity is related to the strong coupling scale of the theory. A simple way to see this 
is by couphng the gauge field to 4D matter sectors localized on the branes, say, to massless 
fermions. At tree level, i-channel amplitudes for two-by-two scatterings of the brane fields 
depend on the scattering energy as M.uv,iR ~ E'^Puv,ir{~E'^)- Below the KK scale the 
theory is efi'cctivcly four-dimensional, thus the amplitudes do not grow with energy. One 
finds M.\]v,iR ^ 9I/ L = Qq and, of course, we assume that the zero mode gauge coupling qq 
is perturbative. Once we cross the KK scale the IR brane-to-brane propagator switches to 
a 5D behaviour. The IR amplitudes then grow as M.ir ~ aJ^Egl and violate the unitarity 
bound at ~ ^s{L) = ai/gi- Above As{L) the IR brane fields are strongly coupled and 
the perturbative description of the IR physics is no longer valid. On the other hand, for 
Mkk < E < k the UV amplitudes evolve only logarithmically, M.\jv ~ glk/ \og{k/ E) (we 
assume that the amplitude remains perturbative in this energy regime) . The linear growth 
starts only above the curvature scale: ior E > k we find Adw ~ Eg'^ leading to the strong 
couphng at -E ~ A5(0) = l/gfy 

We can generalize the preceding arguments by inserting 4D test matter fields at an 
arbitrary position in the fifth dimension. The obvious outcome is that the strong coupling 
scale depends on the position as 

A5(X5) = ^ . (26) 

y5 

The strong coupling scale sets a limit on the validity range of the gauge theory. We are 
forced to cut off perturbative computations below A5 and assume that some UV completion 
(or a non-perturbative formulation) properly describes the physics above A5. Since the 
maximum validity range depends on the position, it is natural to consider a position- 
dependent cutoff, A(a;5) = 0(2:5) A, with A < l/s'i- Deconstruction provides a framework 
to introduce such a cutoff in a gauge-invariant way. 

We have argued that UV brane physics remains perturbative at energies much higher 
than the strong couphng scale on the IR brane. This is possible because of the locahty 
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of higher dimensional gauge theories. At virtuahties E bigger than the KK mass gap, 
the gauge bosons propagate only a distance of order into the bulk, and simply do 
not reach the strongly coupled region close to the IR brane, as can be seen explicitly in 
the approximate Euclidean brane-to-brane propagator of (22). In fact, for E > Mkk the 
Euclidean UV brane propagators would remain essentially unchanged if the IR brane were 
removed, that is for L — > oo. 

There is however a caveat. For time-like momenta, the UV brane propagators have 
poles at the position of the KK masses. At the quantum level, these poles will correspond 
to resonances with finite widths. Measuring these widths, a UV observer can determine 
whether there is strong coupling somewhere else, in particular in the IR brane. This requires 
at least an energy resolution of the order of the KK mass splittings. We can understand 
this in an alternative way by looking at the propagator between the two branes, which is 
no longer exponentially suppressed for time-like momenta, even in the regime p ^ kuL. 
Instead, it has an oscillatory behaviour, so that in principle a UV brane observer is able 
to probe the strong coupling of the IR brane physics. In practice, however, because the 
UV scattered particles have non-zero energy spread F, one must average over this range of 
energies, and the interference of the different modes leads to an exponential suppression. 
This effect can be taken into account by evaluating the propagators at = 2E + iT, 
which leads to a suppression factor ~ exp(— F/oifc) [19]. Hence, we see again that the 
strongly coupled IR brane physics remains screened from the UV brane observers, as long 
as energy resolution is worse than the spacing between KK modes ~ ol/c. 



3 Deconstruction 

We move to the deconstruction of warped gauge theories [9-11]. Our goal for this section 
is to perform an analytical computation of the spectrum and propagators in a model 
approximating the gauge theory in the Randall-Sundrum background. A detailed physical 
discussion of the results is postponed to Section 4. For definiteness, we restrict in the 
following to deconstructing 5D theories with U{1) gauge group. For U{n) or SU{n) the 
analysis is very similar and poses no additional technical problems. 

We consider a U (1)^+-'^ gauge theory with N complex scalars (the links) with charges 
(lj_i, — Ij). The action is given by 

S.^ j d'x ^^^^ + E + - 4)^.f + ^(l^^f)) • (27) 

The gauge coupling of each gauge group is set equal to g. Moreover, at tree level we forbid 
kinetic operators F^^F^^^ with j ^ k or higher- dimensional operators coupling different 
links. These arbitrary choices are meant to reproduce certain features of the 5D theory 
such as 5D coordinate invariance and locality. Furthermore, we assume that the scalar 
potential is such that all the links acquire vev's, {^j) = Vj/\^. One real component of 
each $j gets mass of order vj while the other stays massless (it is the Goldstone boson 
eaten by the gauge field that becomes massive). We isolate the massless components and 
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ignore the massive ones by going to a non-linear parametrization, $j ^^e**^^/^^ , in which 
the action depends only on derivatives of Gj. The action for Aj^ and Gj becomes 




We now compare this action to that of a latticizcd 5D warped gauge theory^. We 
divide the fifth dimension into intervals of size A = (called the lattice spacing) and 
thus introduce N + 1 lattice points yj {yo corresponds to X5 = 0, while yN corresponds to 
X5 — L). We thus rewrite the continuum action (3) as 

(29) 

It is natural to identify the inverse lattice spacing with the cutoff scale A of the continuum 
theory. Comparing the 5D latticized action eq. (29) with that of deconstruction, eq. (28), 
we are able to set up the dictionary between 5D warped gauge theories and deconstructed 
models: 





^ J, 








^ gGj 


a{yj)A 


^ gvj, 




^ 9- 



(30) 

Of course this dictionary is ambiguous at higher orders in the lattice spacing A, as the 
latticization procedure is not uniquely defined (for example, discretizing 85 is ambiguous). 
Specializing to the AdSs background we choose the link vev's as 

Vj = vq^ , (31) 

with < 1 for definiteness.^ This leads to the following AdSs-deconstruction dictionary 



k 
A 



log? ^ 



a{yj) ^ =^ aL'r^ q'^ , 

A ^ gv, 

gl ^ 9 /v. (32) 



^We choose an equally-spaced discretization of the coordinate x^. It is possible to choose a different 
latticization, for instance by uniformly discretizing in the conformal coordinates. This would match a 
different, non-equivalent deconstruction model with the same low-energy limit. 

^See ref. [13] for a discussion of the circumstances under which such a pattern can arise dynamically. 
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The UV brane corresponds to the site j = 0, and the IR one to j = N. UV and IR brane 
matter fields can be represented in deconstruction as matter fields transforming under the 
0-th or A^-th group, respectively. 

One often considers a formal limit of sending the lattice spacing to zero, that is A ^ oo. 
This is called the continuum limit. In deconstruction, this limit corresponds to g — > 1, 
N ^ oo, V ^ oo with and gvlogq~^ kept fixed. Strictly speaking, perturbative 
deconstruction can never reach the continuum limit, which is just a rephrasing of the 
triviahty problem of 5D gauge theories. Indeed, from the last relation in eq. (32) this 
would require either — > oo or — > 0. For this reason we will refer to the continuum 
limit in the following sense: on the 5D side, the continuum theory means the 5D theory 
with large enough A to neglect cutoff effects, in particular A ^ A;, but still A < l/gi- On 
the deconstruction side the continuum limit amounts to choosing 1 — q ^ 1, A^(l — q)>l. 

Previous analytical studies in deconstruction were restricted to the case g <S 1, far 
away from the continuum limit. As can be seen from eq. (32), this corresponds to the 5D 
theory with the cutoff smaller than the curvature scale. Such field theory is UV-sensitive as 
higher dimensional operators in the 5D action, e.g. -^F^jiy, may be sizable. Thus, although 
deconstructed models with q <^ 1 are perfectly well-defined and interesting in their own 
right, their relation with 5D physics is obscure. In this paper we extend analytical studies 
to arbitrary q, including q ^ 1. We will thus be able to approach deconstruction of 5D 
theories with k <^ A, which are well under control on the 5D side. 



3.1 Kaluza— Klein approach 

We now turn to computing the spectrum of the theory. The gauge boson mass terms are 
the following 

>C.a.. = E^A^^(A^-A^-^)^ (33) 

We perform the rotation Aj^ — fj^nA^, which brings the mass terms to the diagonal form 
J^raass = '^n=o ^''^ni^^Y ■ The Coefficients fj^n should be viewed as a discretized KK profile 
fniUj)- They satisfy the following difference equations (we define x„ = rrin/gv) 

(q + q~^ - q~^(xnq~^f)fj,n - qfj+i,n - <f^fj-i,n = (34) 
subject to the boundary conditions 

/o,n — f-l,n -I fN,n — /iV+l,n ■ (35) 

The normalization condition reads g~'^Y^^=o f'j,n — ^- Using the dictionary, we can show 
that the difference equation (34) translates to the continuum equation for the KK profile, 
cf. eq. (7). The boundary conditions are obviously the discretized version of the Neumann 
boundary conditions, cf. eq. (8). 

It is easy to find the zero-mode solution to eqs. (34) and (35): 

= = . (36) 
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The wave function of the massive modes can also be found analytically. We first define the 
variable t\j] = Xnq~^ and the function F{t[j]) = fj,n- In terms of these variables eq. (34) 
becomes a q-difference equation 

{q + q-^ - q-h^)F{t) - F{tq-^) - F{tq) = . (37) 

This equation has been extensively studied in the mathematical literature [20-24] and its 
solutions are called q-Bessel functions. Since, perhaps, readers are not so well acquainted 
with these functions, we summarize all their relevant properties in Appendix B. The q- 
Bessel function Jy{t\q^) is defined by the series (B.2). The other independent solution 
Yj,{t]q^) is called the q- Neumann function and is defined in eq. (B.5). The q-Bessel func- 
tions generalize the ordinary (henceforth referred to as continuum) Bessel functions and 
enjoy similar properties (such as integral representations, recurrence relations). In fact, for 
g — > 1 they are simply related to the continuum ones, see eq. (B.8). 

Equation (37) is a special case of the Hahn-Exton equation (B.l) with v = \ and the 
solution is F[z) = AJiit] q^) + BYi{t] q^). Therefore the deconstructed KK profile can be 
written as^ 

fj,n = KoT^ \Xo{.Xn, q'^)Ji{xnq~^; q^) - Jo{xn] q'^)yi{xnq~^; ?^)] ■ (38) 

The ratio of the two integration constants has been chosen such that the first of the 
boundary conditions in eq. (35) is satisfied (the recursion relation (B.ll) is useful to prove 
this). The other boundary condition determines the KK mass spectrum. We find the 
spectrum is given by solutions of the equation 

Mxn] q^)Yo{q-''-'xn; q^) - Yo{xn, ?') Jo(g-^-'a;„; g') = . (39) 

Note the tantalizing formal similarity of eqs. (38), (39) to the continuum KK profiles (9) 
and the quantization condition (11). In the next section we will specify the parameter 
range in which the continuum and deconstructed physics indeed match. 

3.2 Position-space approach 

Much as the continuum theory, deconstruction admits a position-space picture. The decon- 
structed position propagator for the gauge field will be denoted by -P^^(p^). The indices 
J, /c are obvious analogues of the position variables x^,yr, in the continuum propagator. 
Choosing the gauge-fixing term as 

1 ^ 

^ = E (^M^M + ^aMq'G, - g^+^G,+i))' (40) 

^ ^ i=o 

^We could exactly calculate the normalization of the KK profiles using results about q-integrals [22] , 
but here we simply note that the correct normalizations are automatically included in the propagators in 
the next subsection, and can be obtained from them. 
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removes the mixing between Gj and (in the above Gq = = Gjv is understood). In 
this gauge the propagator satisfies D^J^P^^ = i^^^T]^^, where is the kinetic operator 
(Fourier-transformed to 4D momentum space): 

g^D% = {-p^7]^y + (1 - l/OP^^Pu) 5jk + g^v\^^ {Sjk + q^Sjk - Sj,k+i - q^Sj^k-i) V^^u ■ (41) 
The propagator is of the form P^^(p2) = {r}^^ - E^)p. ,^[p2) + P^p._^(p2/^)^ ^^^^^ 
satisfies 

(_p2 + g^yyj^i + q^)) p._, _ g\Y^{p._,^, + g^P.+i,,) = ig%,k , (42) 
together with the boundary conditions 

Po,k = P-l,k , PN,k = PN+l,k ■ (43) 

These boundary conditions, technically speaking, are chosen such that eq. (42) is correct 
for j — and j — N. They are analogues of the Neumann boundary conditions for 
the continuum propagator. See Appendix A for the deconstruction of general boundary 
conditions. 

For j 7^ k, the propagator equation (42) is the same as the KK mode equation (34) 
with Xn ^ X = p/gv. Therefore we can easily write down the solutions for j < k and j > k 
that satisfy the relevant boundary condition: 

P<, - A<q-' [Yo{x; q')Mq-'x; q') - Mx; q')Y^{q-^x; q')] 

P>, = A>q-' [Yo{q-''-'x; q')Mq-'x; q') - Mq-'^-'x; q')Yi{q-^x; q')] . (44) 
The constants A^'^ are determined by the matching conditions 

Pk,k = Pk,k ' Pk+l,k = Pk+l,k + ^2q2k+2 ' f^"^^) 

which make the solution (44) satisfy the propagator equation for j — k and j — k + 1. 
After some algebra we solve for A^'^ and derive the deconstructed position propagator: 

PAp') = ^2^i_q2f ''~' [Yo{x;q')Mq-^x;q') - Mx; q')Y^{q-^ x; q')] 
[Yo{q-''-'x; q^)J,{q-''x; q^) - Mq-^'-'x; q^)Y^{q-^x- g^)] (^6) 

for j < k. The denominator reproduces the poles at equal to the KK masses given 
byeq. (39). 

We can also define "branc-to-branc" propagators 

tg Y,{q-^-^x- q^)Mx; q") - Mq-^'-'x; q')Y,{x- q") 



Pooip') = 
PnAp'') = 



vp Yo{x; q'^)Jo{q-^-^x; q^) - Jo{x; q'^)Yo{q^'^-^x; q^) ' 
ig Yq{x] q^)Ji{q-^x\ q^) - Jo{x; q^)Yi{q-^x; q^) 
qNyp Yo(x; q'^)Jo{q~-'^~^x; q^) — Jq{x; q'^)Yo{q^^^^x; q^) 



p / 2^ ^ - q^) I /.^x 

V Yo{x-q^)Jo{q-^-^x;q^)-Jo{x-q^)Yo{q-^-^x-q^)- ^ > 
13 



We will also need the UV propagator evaluated at Euclidean momenta = —p%. It is 
convenient to rewrite it in terms of the modified q-Bessel functions Iu{x; q^) and K^{x] q^) 
defined in eqs. (B.17) and (B.18): 

ig Ka{q-^-^XE]q^)h{xE]q^) + h{q~^'^XE]q^)Ki{xE]q^) , . 

where xe = PeIqv. 

The connection between the KK and position pictures in deconstruction is provided by 
the spectral formula 

fj,nfk,r, 



PAp')--iY.^¥^^ (49) 

71=0 ^ " 

which involves a finite sum. This implies that the propagator (46) is, in fact, given by a 
ratio of two polynomials in . 



3.3 Holographic approach 

Finally, it is possible to implement the holographic approach in deconstruction. In fact, at 
the conceptual level, constructing the boundary effective action is even clearer here. The 
procedure in deconstruction amounts to integrating out N gauge bosons A^^ with j > 1 
and leaving as the low energy variable.^ We obtain, for ^ = 1, 

1 f 



4 Discussion 

In this section we give a detailed discussion of the spectrum, the KK profiles and the 
propagators derived in Section 3. We first consider the case q ^ 1 and determine the 
energy range in which deconstruction is a good approximation of the continuum 5D gauge 
theory. Then, we extend our discussion to general q, including the case g <S 1, which does 
not have a 5D interpretation. Finally, we comment on the holographic interpretation of 
the deconstructed setup, for both the continuum case g ~ 1 and for smaller values of q. 

Let g ~ 1, but still q^ ^ q^^^ <^ 1. In this case it is convenient to rewrite q = 1 — 5 
with 5 (but not NS) being a small parameter. From eq. (32) 5 is interpreted as the ratio 
k/K, that is k 5gv u-p to corrections of higher order in 5. From this it follows that 
the smallness of 5 is a necessary condition for the corresponding 5D theory to be in a 
controllable regime. 

We first discuss the KK spectrum of deconstruction for 5 <^ 1. For small enough 5, 
according to eq. (B.8), q-Bessel functions can be approximated by the continuum ones, as 

^Note that neither A° nor the remaining gauge bosons are mass eigenstates. However, integrating out 
linear combinations of Ught and massive fields is equivalent to integrating out heavy mass eigenstates after 
appropriate redefinitions of the light fields in the low-energy theory [18]. 
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long as their argument is less than 1. Thus, for m„ < q'^gv, the eigenvalue equation (39) 
can be approximated as 

The spectrum obtained by solving eq. (51) matches the continuum spectrum with k = Sgv 
and gl = , as prescribed by the dictionary (32). In particular, it is linearly spaced: 

rrin ^ 'K5gvq^{n - ^) , "^n ~ q^gv ■ (52) 

One can also easily see that the deconstructed position propagators match the continuum 
ones (more precisely, P{p^ ,y^) ~ Pjkip^) after relating the parameters as in eq. (32)) 
for momenta smaller than gv. This ensures that both theories describe the same physics 
below the scale q^^^gv at leading order in 5 and p/gv. 

For rrin > q^gv, some of the q-Bessel functions describing the KK profile and the mass 
spectrum fall outside the range of continuum approximation. In particular, the functions 
YQ^rrin/q'^'^^gv; q^) and jQ{mn/q^~^^gv; (f) in the quantization condition (39) can no longer 
be approximated by continuum Bessel functions. According to eq. (B.34) these functions 
oscillate rapidly with exponentially decaying/growing amplitudes. On the other hand, 
for m„ <^ Sgv, we have \YQ{mn/ gv\q^)\ ^ \Jo{jiin/ gv\q^)\. Wc can thus infer that the 
solutions to eq. (39) for q^ gv ^ m„ < 5gv are approximately given by the asymptotic zeros 
of Jo{mn/q'^~^^gv; (f). Using, eq. (B.35) we find an exponential spectrum, rrin ~ gvq^~^^q~"', 
in this regime. This formula can be refined using the detailed asymptotic behaviour of the 
q-Bessel functions at small and large argument, given by Eqs. (B.26), (B.34) and (B.41). 
We obtain^ 

rrin ~ gvq^^^q'"-^"'" , q^^^gv < rUn ^ Sgv , (53) 

where 

The q-Euler-Mascheroni constant 7^2 is defined in (B.27). The upper limit in the region 
of validity of (53) arises from the use of the leading terms in the expansions of the q-Bessel 
functions at small argument. Moreover, in order to be able to use (B.41), we have made 
the ansatz |q;„| ^ 1, which is indeed satisfied by (54) when m„ < Sgv. 

Closer to the top of the KK tower, for Sgv < m„ < 2gv, we have not found an accurate 
analytical formula for the spectrum. But from the asymptotic form of the q-Bessel functions 
and the trend for the masses below Sgv, it can be argued that a growth stronger than 
exponential will persist until the scale ~ ^gv. This agrees with an extrapolation from 
the known results at small q and with the exact formula for the product of non-vanishing 
eigenmasses in Ref. [12]. The numerical computations confirm this behaviour. For m„ > 
2gv, there is no solution to the quantization condition: the deconstructed KK tower is cut 
off. 

''Note that for weak warping, < i5, there is no regime where this approximation is valid. 
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Figure 1 : Numerical comparison of the deconstructed KK spectrum (red circles) and the continuum one 
(black line). KK masses are given in units of k. The deconstructed spectrum is calculated for N = 75 and 
q = 0.97, which corresponds to the warp factor ql ~ 0.102 on the continuum side. 



Summarizing, we see that above the scale q^gv, the spectrum depends (approximately) 
exponentially on n, rather than linearly, and resembles nothing one encounters in a 5D 
theory.® It is worthwhile to emphasize that the linear "continuum" regime is not just the 
linear approximation of the subsequent exponential regime, at least for small 6. In fact, 
for the first few modes in the exponential regime, the masses grow more slowly than in the 
linear regime. This can be clearly seen in Fig. 1. 

A similar analysis can be applied to study the KK profiles. As we have already men- 
tioned, for KK modes with masses below gvq^ they match the continuum ones. Above 
this scale, the deconstructed KK profiles are still described by continuum Bessel functions 
for j < j*, where q^' = rrin/gv. However, since the mass of the deconstructed n-th mode 
is different from the continuum one, the deconstructed profile is different as well. The 
normalization factor is different too, because of the behaviour at large j. For j > j^,, 
the quantization condition (39) and the asymptotic behaviour of the q-Bessel functions 
given by eqs. (B.32), (B.34) and (B.41) imply that the deconstructed KK profiles are ex- 
ponentially damped, in contrast with the continuum ones, which oscillate with increasing 
amplitude. Hence, above gvq^ , deconstructed KK modes are decoupled from the IR brane, 
while they are more strongly coupled than the continuum ones to the UV brane, because 
of a larger normalization factor. We plot in Fig. 2 a typical deconstructed KK profile in 

^The exponential spectrum (53) is a specific feature of the deconstructed model considered here. We 
have checked that the spectrum outside the continuum regime is very different in the deconstructed model 
corresponding to a uniform discretization in the conformal coordinates. 
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Figure 2: Deconstructed (red continuous line) and continuum (blue dashed line) KK profiles of the mode 
n = 55, for N — 100 and q — 0.955, which corresponds to the warp factor « 0.010. 



the regime m„ ^ gvq^ , together with its continuum counterpart. 
The preceding discussion identified the deviation scale 

AD = 9vq'', (55) 

above which the deconstructed KK spectrum and profiles no longer match those of the 
continuum theory. Using the dictionary (32) the deviation scale is translated to Aa^ = 
A(L), which is the cutoff scale on the IR brane. Note that for g ~ 1 this scale is much 
larger than the KK scale Mkk ~ 5gvq^ and there are ~ 1/5 KK excitations that are 
properly matched. 

However the applicability range of the continuum theory is far larger, as long as we 
keep to the appropriate observables. For example, as we discussed above, 5D UV brane 
physics remains perturbative up to a much higher scale, of order A5(0) = 1/(75. Naively, it 
would seem that deconstruction cannot reproduce continuum UV brane physics above the 
scale Kd- We will show however that it does. 

Consider the deconstructed UV propagator at Euclidean momenta, eq. (48). According 
to (B.36) and (B.45), the modified q-Bessel functions eq. (48) share the property with their 
continuum analogues that at large argument K^{x;q^) <^ Iy{x;q^), including the 1 
region where they deviate from the continuum K and I functions. Thus, for pe ^ 6gvq^ 
we can approximate 

vpEKo{^;q^) vpE Ko{§^) 
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Figure 3: Deconstructed (red continuous line) and continuum (black dashed line) Euclidean brane-to- 
brane propagators at the UV (left) and IR (right) branes, for iV = 70 and q = 0.95, and thus ~ 0.0276. 
Momenta are given in units of k. 



The last equality involving only continuum Bessel functions holds for pe ^ gv. We can see 
that Poo correctly reproduces the continuum UV propagator (with k = 6gv, gl = g/v) at 
energies below gv, cf. (21) and (22). This ensures that the two frameworks describe the 
same UV brane physics (for example, scattering amplitudes of UV brane localized matter 
fields) below gv. The limiting scale is translated to A = A(0) - the cutoff scale on the UV 
brane. 

On the other hand, on the IR brane the matching is terminated at a lower scale. For 
Pe ~ gvq^ the deconstructed propagator is purely four-dimensional, PNNi—PE) ~ w"^ Ip\i 
and deviates from the 1/pe behaviour of the continuum IR propagator. Thus IR physics 
is matched by deconstruction only below the scale gvq^ , which is translated to the cutoff 
scale on the IR brane A(L). We compare the continuum and deconstructed UV and IR 
propagators in Fig. 3. 

This discussion can be easily generalized to show that, for an arbitrary position y in 
the bulk, deconstruction reproduces continuum physics up to A{y) = Aa{y). Therefore, 
the deconstruction framework we consider is indeed a realization of a position-dependent 
cutoff. 

Let us pause for a moment to summarize what we have shown so far. Our results imply 
that deconstruction provides a correct approximation of the continuum physics at the UV 
and IR branes, and at any position in the bulk. The matching holds all the way up to 
the respective position-dependent cutoff scale. In particular, for A ~ l/^f^ deconstruction 
works well all the way up to the position-dependent strong-coupling scale, that is in the 
entire perturbativity range of the continuum theory.^ The correspondence extends to the 
energy range where the KK spectra and profiles in the two theories are completely different. 

^On the other hand, if we choose A sufficiently smaller than the 4D behaviour of the deconstructed 
theory sets in early enough to stop the linear growth of tree-level amplitudes and prevent the occurrence 
of strong coupling. Therefore, this scenario can provide, in principle, a perturbative completion of the 
continuum AdSs theory that can be valid up to much higher scales. 
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We were able to establish the agreement of the deconstructed and continuum theories 
in the position approach. However, it all seems miraculous from the point of view of the 
KK approach. As we have discussed before, the deconstructed KK masses and profiles 
deviate from the continuum one at the scale A/). Thus, in the two frameworks scattering 
amplitudes above Ad involve an exchange of a different number of KK modes with different 
masses and couplings to the brane fields. Somehow, for UV brane amplitudes the two effects 
cancel and both frameworks predict the same result. The reason for the better agreement 
of the UV propagators is the locality of the action (28) in the extra-dimensional lattice: 
the Euclidean UV propagator is screened from the IR brane. On the other hand, the KK 
masses are determined by global properties of the theory and are thus sensitive to the 
physics at the IR brane. 

Of course, strictly speaking, the matching of the propagators holds only for t-channel 
amplitudes, when the propagators carry space-like momenta. An s-channel experiment on 
the UV brane with enough precision could resolve individual KK resonances, which do not 
fit in the two theories, and distinguish our completion from a strongly-coupled continuum 
theory. In this case, above Ad, the matching holds only when the energy spread of the 
scattered beams is sufficiently larger than the KK spacing, so that individual resonances 
cannot be resolved. 

We now move to discussing the situation when the parameter q is not close to 1, 
that is when deconstruction is away from the continuum limit. In this case there is no 
energy regime where the deconstructed spectrum mimics the continuum one. However, 
earlier works that concentrated on g <^ 1 showed that also in this regime deconstruction 
reproduces certain features of the continuum theory, for example logarithmic running of 
gauge couplings. We are now in a position to make the relation more precise and interpolate 
between the different values of q. 

At the level of the spectrum, for general q we can again distinguish three regimes. For 
TTT-n ~ q^^^gv, the spectrum is also linear, but now this regime holds only for a few modes. 
For nin ~ q^^^gv, the KK masses are given with great precision by (53), and only the last 
few modes with (1 — q)gv < m„ < gv deviate; see Fig. 4. In the limit g <i; 1 we see that 
only the exponential regime survives, in agreement with [10-12]. The behaviour of the KK 
profiles is qualitatively the same as the one discussed before. 

Consider now, once more, the UV propagator (48), but now for arbitrary q. At 
large argument / still dominates over K, so that at energies above the KK mass gap, 

~ gvq^'^^i we can approximate 

vpEKQ{pE/gv]q^) 

As long as pe ~ gv, we can use the small argument asymptotics (B.26) of the modified 
q-Bessel functions to obtain 

p (_^2 ^ ^ 1 ig'^^ogq-^ ...^ 

-fool Pe) ~ :j — 7 — i\~=777T75 — V ■ \^°) 

log {gv{l - q^)e ^"^qy^PE) Pe 

The deconstructed UV propagator exhibits analogous momentum dependence to that of 
the continuum UV propagator below the curvature scale, cf. eq. (21). The factor g^k in the 
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Figure 4: Comparison of our approximate analytical spectrum (black empty diamonds) and the numerical 
one (red solid circles) for = 50 and q = 0.832, and thus « 10~*, in units k = 1. We plot only the 
first 40 modes, which lie in the region of validty of our approximation. The relative error is below 1% for 
all but the last few modes. 



continuum is replaced here by log q~^, in agreement with the dictionary (32). Away from 
the continuum limit, for q much smaller than 1, only the argument of the "classical log" 
does not match the continuum value. Therefore deconstruction with arbitrary q captures 
the essentials of the continuum UV brane physics below the scale k. This explains why 
deconstruction could reproduce certain features of the 5D theory, such as the logarithmic 
running of gauge couplings [25], also in the parameter space away from the continuum 
limit. 

This brings us to a comment on the relation of deconstruction and holographic CFTs. 
We have argued that the UV boundary physics is reproduced in deconstruction for arbitrary 
q, below the scale gv where the approximation (58) is valid. Therefore it is natural to 
conjecture that deconstructed models, also in the parameters range where they do not have 
a 5D interpretation, describe the dynamics of some large Nc strongly-coupled theories which 
are approximately conformal over a range of scales. This is an interesting generalization of 
the AdS/CFT conjecture. In the case at hand 4D strongly-coupled theories might be dual 
to weakly-coupled 4D theories with extended gauge symmetry. Matching the propagator 
(58) onto the quark bubble calculation in QCD at large Euclidean momentum, the number 
of colours in the CFT can be related to the parameters of deconstruction: 



127r2 



(59) 
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^ log q 



1 ■ 
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According to the dictionary (32) we have g^k < — > g^logq ^, so this a translation of the 
well-known AdS/CFT relation A^^ = Un^/glk. 

5 Summary and outlook 

In this paper we clarified in what sense and in which parameters range deconstruction 
approximates 5D gauge theories in the Randall-Sundrum background. In our analysis 
of the deconstructed theory, we employed powerful tools of the mathematical theory of 
q-Bessel functions. This allowed us to study analytically the parameter space that was 
previously accessible only to numerical methods, including the ~ 1 region where the 
corresponding 5D theory is under control. 

The main result of this paper is an explicit proof that a warped 5D gauge theory 
can be approximated by deconstruction in all its perturbativity range. More precisely, 
deconstruction provides a faithful description all the way up to the position-dependent 
cutoff K{x^) = 0(2:5) A, where A is the inverse lattice spacing. In particular, the continuum 
theory with A ~ l/g^l can be reproduced by deconstructed models with perturbative gauge 
coupling, ~ 1, in which case the matching extends all the way up to the position- 
dependent strong- coupling scale A5(a;5) = a{x^)/g1. This is not a trivial result, as the KK 
spectra of the two theories deviate at a much lower scale, A^ = o^A. 

The technical results we derived can be readily apphed to more phenomenological stud- 
ies. In fact, we have shown that computations in deconstructed warped theories can be 
performed at the quantitative level comparable to that in 5D. This is an encouragement 
to study deconstructed versions of phenomenological models in AdSs, such as the models 
of the electroweak sector of refs. [2]. Deconstruction could provide for a UV completion 
of these models that could remain perturbative to much higher energy scales. Further- 
more, deconstruction offers a playground to study the evolution of gauge couplings in 
AdSs [17,25]. This was already exploited in [10,11], but only in the region g <C 1, where 
the link with the 5D computation is obscure. The virtue of deconstruction in this case it 
that it provides a concrete physical realization of the 5D cutoff physics, which allows, in 
particular, to study threshold effects. 

It would also be interesting to exploit the relation of deconstruction with holographic 
CFTs. Most interestingly, deconstruction can describe certain aspects of the real-world 
QCD in the strongly coupled regime. For example, chiral symmetry breaking and physics 
of p and a resonances can be captured, as shown in ref. [26]. More recently, effective 
description of low energy QCD was pursued in 5D continuum models [7] . Our results allow 
for quantitative studies of deconstructed versions of these models. This could shed further 
light on the origin of this AdS/QCD correspondence and its connection to the Migdal's 
approach [27] discussed recently in [28]. Indeed, the high energy behaviour of the continuum 
correlation functions, which is non-analytic in p^, is approximated in deconstruction by a 
ratio of polynomials, similarly as in the Migdal's approach. 

Finally, we note that the methods used here can be directly applied to fields with 
different spins, at least in the masslcss case. For instance, the difference equation for the 
graviton corresponds to the Hahn-Exton equation (B.l) with v — 2. Therefore, all the 
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results presented here for gauge bosons translate in a straightforward way to the gravity 
case. In particular, the qualitative features of the spectrum are the same, and agree with 
the discussion presented in ref. [15]. It is also hkely that the methods of q-difference 
equations can be applied to study deconstruction of more general backgrounds than the 
Randall-Sundrum AdSs one. 
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Appendix A Implementation of general boundary con- 
ditions 

We generalize our results to gauge theories with Dirichlet boundary conditions on one or 

on both branes. 

The mixed momentum-position space propagator for a gauge boson in a slice of AdS^ 
can be succinctly written as 

(A.l) 

where the parameters A"[0], i^"[0], ^'^[1'], Bf^[L\ depend on the boundary conditions. They 
take the values 

if Neumann boundary conditions are imposed at X5, and 

for Dirichlet boundary conditions at x^. We can also define brane-to-brane propagators, 
as long as we deal with Neumann boundary conditions on the given brane. For example, 
for the Neumann-Dirichlet case we obtain the UV brane-to-brane propagator 

P,y{p)-P (^'O'O)- pyo(f)Ji(^)-Jo(f)n(^)' ^^-'^ 
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which is suppressed for yj—p^ <^ M^^k-, while for yj—p^ ^ M^^k it is indistinguishable 
from the UV propagator in the Neumann-Neumann case. 

We turn to discussing the method of realizing these more general boundary conditions 
in deconstruction. To this end we include two more scalar fields: $o charged under 
with a vev vq/\/2 and ^n+i charged under with a vev fTv+i/y^, so that the boundary 
entries in the gauge boson mass matrix are modified. The effect of the additional vev's is to 
modify the "boundary conditions" for the propagator. More precisely, eq. (42) represents 
the correct propagator equation if we impose 



PN+l,k — PN,k ■ 

Setting t>o = f at+i = we recover the Neumann boundary conditions (43). With vq — vwe 
get = which mimics the Dirichlet boundary conditions on the UV brane. Similarly 
vn+1 — q^~^^v yields PN+i,k = 0, which is a deconstructed version of the Dirichlet boundary 
conditions on the IR brane. Other choices of vq and vn+i correspond to mixed boundary 
conditions on the continuum side. 

We are now in a position to write the general form of the deconstructed position prop- 
agator: 

P<^P(j,2\ ^J, -k-j 

[A-[0]Mq-^x; q^) - B'-[Q\Y^{q-^ x- q^)\[A^[L\J^{q-^x- q^) - BP[L\Y^{q-^x- q^)\ 

A<^[Q\BP[L\ - B<^[Q\A^[L\ 



(A.6) 

Neumann (N) or Dirichlet (D) boundary conditions specify the parameters in eq. (A.6) as 
follows 

A''[Q]=Yo{x-q^) i?^[0] = Jo(x;g^) 

/1^[L] = Y,{q-^-'x- q') i?^[L] = Mq-''-\ q') 

A^[0] = Yi{qx; q^) S°[0] = Ji{qx; q^) 

A^'lL] = Yi(?-^-ix; g2) ^^[L] = Ji(g-^-ia;; q^) . (A.7) 

The correspondence between eqs. (A.l) and (A.6) can be established using the methods 
discussed in Section 4. In particular, below the deviation scale A/j = gvq^ (translated 
to the IR cutoff scale A(L) = clA), physical amplitudes in deconstruction mimic those 
in the continuum theories at the leading order in 5 and p/ gv. Above the deviation scale, 
the spectra (the poles of the propagator) are different. However t-channel amplitudes on 
the UV brane are still reproduced in deconstruction, all the way to the UV cutoff A(0) = 
A. Therefore for general boundary conditions the correspondence between deconstruction 
and continuum holds at the same quantitative level as for the Neumann-Neumann case 
discussed in Sections 3 and 4. 
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Appendix B Q-tutorial 



This appendix contains definitions of the q-Bessel functions and a review of their vital 
properties. Our presentation is based mainly on ref. [21]. We also derive several results 
concerning the asymptotic behaviour of the q-Bessel functions, which have not been given 
in the literature. 

Several inequivalent q-analogues of the Bessel functions have been studied in the math- 
ematical literature. The ones relevant to our purpose are solutions of the so-called Hahn- 
Exton q-difference equation, 

{q^l'^ + g-'^/' - q-''/h^)F{t) - F{tq-^'^) - F{tq^'^) = . (B.l) 

For definiteness we consider < g < 1. One solution of this equation is the Hahn-Exton 
q-Bessel function [20] (q-Bessel in the following) denoted by Ji,{t;q). It is defined by the 
power series 

where the q-shifted factorials are defined as 

for a e C and k e Z+ = {0, 1,2,.. .}, and (a; g)oo = linife_*oo(o; Q)k- The q-Bessel Ji,{.; q) is 
analytic in C\{0}. For u — n & 'Z and \t\ < 1, it has an integral representation: 



27r Jo (te*^; q). 

It is known that all the zeros of the q-Bessel function of order i/ > — 1 are real and that 
the non-zero ones are simple [22]. 

For u ^ Z Jv{t]q) and J-i,{tq''^/'^;q) are two independent solutions of eq. (B.l); how- 
ever, for integer u = n, there is the relation J^n{t;q) = (— l)"g2 J„(tg2"; g). Therefore 
one introduces^° the q- Neumann function Yi,{t;q), which is an independent solution for 
arbitrary v. 

Y,{t-q) ^ EiMEi(i^^-4 [cos(7rz/)gt J,(t;g) - J_,(tg-t;g)] , (B.5) 



TT 

where 



^°Our definition differs by a factor q "2 from the one in [23]. With our definition, J and Y satisfy the 
same recurrence relations. 
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is a q-extension of the Euler F- function, satisfying limq_^i- Tg^u) = r(i/). For integer u, 
the hmit — > n e Z is understood on the r.h.s. of eq. (B.5). Exphcitly, this hmit gives [23] 



TT log q 1 — 9 



n-1 , ^ .2fc 



1 - ^ 9)n-fc-lt' 



^log? f::t (?;?)fe(?;g)n+fe lrg(n + A; + 1) rq{k + i) 



k=0 

oo 



^ {q;q)k{q;q)n+k 

for n G Z+ (omitting the second term when n = 0). For negative integers Y'_„(t;g) = 
(-irgtF„(tgt;g). 

We now review properties of the q-Bessel and q-Neumann functions. Since the form 
relevant to deconstruction is J„{t;q'^) and Yi,{t;q'^) (rather than Ji,{t;q) and Yi,{t;q)), we 
find it convenient to present all formulas in this form. 

The q-Bessels and q-Neumanns are q-analogues of the ordinary Bessel and Neumann 
functions, in the sense that there exists the "continuum limit" g — > 1, 

lim Z,{{l-q)t;q'')^Z,{t), (B.8) 

with Z = J,Y and J,y{.),Yi,{.) the ordinary Bessel and Neumann functions. For q close to 
1, the continuum approximation holds for (1 — -C 1, as can be seen from both the 
series and the integral representations. In fact, when 5 = 1 — 5 ~ 1 and \t\ <^ 1, we can 
Taylor expand Flq'^H) in powers of {5t). Keeping terms up to the second derivative, we 
obtain 

SH^F"{t) + 5HF\t) + {q-^e - - g-")F(i) = , (B.9) 
which is solved by the Bessel or Neumann functions: 

F{t) = Z4sinh^[logW2] « (0 ■ (B-10) 

This derivation is valid as long as higher derivatives in the Taylor expansion of F{q^^t) can 
be neglected, that is when \6tF"'{t)\ <^ \F"{t)\. For the solution eq. (B.IO) this is indeed 
true for |i| <^ 1. 

The q-Bessels and q-Neumann satisfy recursion relations analogous to the ones of their 
continuous cousins. In particular, the following ones will prove very useful 

q''Z^{t; q^) - Z^{qt; q^) = -qtZ^+i{qt; q^) , 

Z,{t- q^) - q'Z^iqt- q^) = tZ,-i{t- q^) , (B.ll) 
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with Z = J,Y. These are actually discrete versions of difference-recurrence relations. 
Indeed, one can define the q-derivative of a function as 

DJ(t) = l '"^^ , (B.12) 

which is a q-analogue of the continuum derivative: liniq^i- Dqf{t) = fit). Then, (B.ll) 
are a particular case {k — 1) of [21]: 

[t-'D,] ' [t^Z,{t; q^)] = q^) . (B.13) 

For q ^ 1 this reduces to the familiar recursion relations of the continuum Bessel functions. 

The Wronskian of a differential equation has also its q-analogue. In general for a q- 
difference equation of the form F{qt) + F{q~H) + P(t, q)F{t) = 0, its two independent 
Fi{t) and F2{t) solutions satisfy the relation Fi{qt)F2{t) - Fi{t)F2{qt) = C{q) with the 
r.h.s. independent of t. For q-Bessels, the q- Wronskian relation is 

Mqt; q')Y,{t- q') - Y^{qt; q')Mt; q') = ^"^^^ ~ , (B.14) 

TT 

which can be proved using the scries expansion of the q-Bcsscls. The fact that the q- 
Wronskian docs not vanish at any point t implies that a general solution of the q-difference 
equation (B.l) can be written as [23] 

F{t) = A{t)Mt; q') + B{t)Y,{t; q^) , (B.15) 

where A{t) and B{t) are q-periodic, i.e. they satisfy A{qt) — A{t), B{qt) — B{t). But since 
wc arc interested only in discrete values of the argument, to, qto, qHo, . . ., we can simply 
treat A and B as constants. 

We can also define q-analogues of the modified Bessel and Bessel-Macdonald functions, 
which are solutions of the Euclidean version of the q-difference equation (B.l), 

+ ^ q-u/2^2^p^^^ _ F(^tq-i/2^ _ F{tq^'^) = . (B.16) 
We define the q-modified-Bessel as^^ 

h{t;q)^i-^Mit;q), (B.17) 

and the q-Bessel-Macdonald as 

K^t; q) = ^g-^r,(i/)r,(l - u) [l-.{q-h; q) - q^^t; q)] , (B.18) 
^^Here and in the following i represents e'^'/^ when — tt < argf < 7r/2. 
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with a limit understood for integer v. Our definitions agree with the ones in [24], up to a 
change of variables and a global factor. Equivalently, we can write 



TT 



TT 



(B.19) 



Taking the limit of integer index, we find 



TT 



for n e N and 



logg(gi ";g)n-i 
l-q) 



TT 



Ko{t;q) = -i 



-logg 



Jo{it;q) +iYQ{it;q) 



The continuum limit is 



lim X,{{l-q)t;q^)^X,{t), 



(B.20) 



(B.21) 



(B.22) 



with X = I,K. The q-modified-Bessels and the q-Bessel-Macdonalds satisfy the recurrence 
relations 



[r'D,]'[t-''X^(t;q')] 



X,+k{tq';q') 



[t-'D,]'[t''X^(t;q')]^(-l) 



ak 



4-1'— k 



:Xi,_k(t] q^) , 



(B.23) 



with q; = 0, 1 for X = 7, K, respectively. The q-Wronskian is 

h{qt; Q')KAt; q') - K^iqt; q^)!.^ q') = -^lllz^ 
and a general solution of (B.16) can be written as 

F{t) = A{t)h{t; q^) + B{t)K,{t; q^) , 



(B.24) 



(B.25) 



with A, B q-periodic functions, which again we can consider constant. 

The behaviour of all these functions at small values of the argument can be read from 
the corresponding series expansions. We list below the small-argument asymptotics that 
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are useful for our computations 

Jo{t;q^)^Io{t;q^)^l, 



Yo{t; q^) ^ log + 7,2 

— TT log g \ 1 — q^ 

Y,{t;q')^-^t-\ 

TT 

(B.26) 

for \t\ <^ 1. We have defined the q-Euler-Mascheroni constant 7^ = —1^(1). Its series 
representation is 

00 2k 

7,2 = log(l - q') - log q' ■ (B-27) 

k=l ^ 

Varying q from to 1, 7,2 grows approximately linearly from to 7 = 0.577 .... 

Obtaining the asymptotic behaviour at large argument is more involved. We first look 
at the asymptotic form of the Hahn-Exton equation (B.l). Writing F{t) = Il{t)G{t), 
where 



U{t) = l[{l-q'H')^{qH';q'U (B.28) 

k=l 

and using U{qt) = n(t)/(l - qH'^), U{q-H) = n(t)(l - t^) the Hahn-Exton equation takes 
the form 

G{qt) + (1 - t^){l - qH^)G{q-H) + (1 - qH^){t\-^ - q" - q'^Git) = , (B.29) 

which for |i| 1 reduces to 

G{qt) + qH^G{q-H) - qH^q-^Git) = . (B.30) 

One solution follows from neglecting the first term: G{t) = const t". From eq. (B.30) we 
see that this approximation should work fine very quickly when |t| becomes larger than 
1. To find the second solution it is more convenient to insert F{t) = (l/n(t))G(t) in eq. 
(B.l), which leads to the equation 

G{qt){l - t^){l - gV) + Giq-H) + (1 - f){fq-'' ^ q" - q-^Git) = , (B.31) 

approximated by G{qt) = q~^~^G{t) = at large \t\. The solution is G{t) = const 
Therefore, a general solution to the Hahn-Exton equation has the asymptotic form 

F{t) « At^'iqH'; q'U + B , ' , \t\ » 1 . (B.32) 

[q 1^ iq )oo 
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For imaginary t one solution grows rapidly (faster than any power of t) while the other 
decays (equally rapidly). For real t the situation is more complicated as the "decaying 
part" has poles at tn — q~"'- In the following we pinpoint the coefficient A for the q-Bessel 
functions defined previously. 

We can obtain the asymptotic behaviour of the q-Bessel by approximating (g*^^^; q)k ^ 
{q''~^^',q)oo in the power series (B.2). This is justified because for t 1 the series is 
dominated by terms with large k for which this approximation is fine. Using then the 
so-called q-binomial theorem (see for instance the second reference in [20]), 



n=0 



qn(n—l)^2n 

{q';q% 



we arrive at 



lil > 1. 



(B.33) 



(B.34) 



{q^;q^)oo 

This fixes the normalization of the growing part of the q-Bessel. From this expression we 
can directly read the asymptotic zeros of the q-Bessel J^{t; g^), which are located at 



tr ^q~\ r en, g"'' > 1 



(B.35) 



independently of the index v (though, in fact, the approximation (B.34) is better for larger 
values of u). Since the above derivation is valid for any argument of t, we also find 



h{t- q') ^ t 



{q^;q^)c 



Itl > 1. 



(B.36) 



Finding the asymptotic form of the q-Neumann is even more involved. For integer index 
n e using (B.5) and (B.34) we have 



Ynit;q') 



(-l)"(l-g^)(g^gVi 

27rlogg(g2;g2)^(^2-2n.^2)^_^ ^ 

Yj+21ogg2^ 

r=l 



q^^'-^hog{qt)e{qH^-q\ 



log 



,2(r-n)^2 



t-^{q^-^'H';q^) 



ltl > 1, 

(B.37) 



for neN, 
Yo{t;q'): 



TTlogg ^(g2;g2)^ 



f 1 °° 

logt + logg i^-^Y 

\ r=l 



2J2r 



rq 



(gV;g%, |i|>l 



(B.38) 

Note that the poles in the sum cancel against zeros in (g^~^"i^; g^)oo- For real t, writing 
t — g~*+", with s e N and a e] — 1/2, 1/2], and assuming \a\ <^ 1, we can approximate 
the sum at large t by 



^2g2(r-n) 

2^ I- q^{r-n)l2 - log q 



r=l 



Q~2a _l V / 



(B.39) 
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Therefore, 

(-ir(l-g2)(g2.^2)^^_^ 



27rl0gg(g2;g2)^(^2-2n.^2)^_^ 

1 



1)" ( logi + 2 ( 1 + -3^^ ) logg ) q\ 



, |q;| < 1, i > 1, 
(B.40) 

for n e N, and 



^o(^; q') :^l^7^ + ^=^1 ) (^"t'; q'U , |a| « 1 , t » 1 . (B.41) 

The corresponding expressions for imaginary argument t = ir, r G M, are simpler. In 
this case we can evaluate the relevant sum to leading order at large |r| without any further 
restriction: 

~ T-2„2(r-n) 1 2 i 

-V =!2^-n+- + 0(T-') , (B.42) 

l + g2(r-n)^2 2 log Q 2 ^ ^ ' ^ ' 

r=l 

which leads to the asymptotics 

" ' '^ii^^^SS^i^n' 

for n e N and 

Yo{it;q-')^i (-g'^';g')oo, teK, (B.44) 

2iogg nr;r)oo 

This allows us to build a linear combination of the q-Bessel and q-Neumann, which is purely 
decaying at large imaginary values of the argument. In fact, from the above expression we 
see that this combination is proportional to the q-Bessel-Macdonald. In other words, 

Kn{t;q)<^l, teR, i>l. (B.45) 

Therefore, the q-Bessel-Macdonald shares with its continuum counterpart the property 
that, at large real argument t, it is damped. This is true in the "continuum" region as well 
as for t ^ 1 where the continuum approximation fails. We also see that in this regime 
Kn{t]q) ^ In{t]q)- It is plausible that the same behaviour holds for non-integer indices, 
though we have not proved this result. 
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